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Abstract 

Within the Dijkgraaf-Vafa correspondence, we study the complete 
factorization of the Seiberg-Witten curve for U{Nc) gauge theory with 
Nf < Nc massive flavors. We obtain explicit expressions, from ran- 
dom matrix theory, for the moduli, parametrizing the curve. These 
moduli characterize the submanifold of the Coulomb branch where all 
monopoles become massless. We find that the matrix model reveals 
some nontrivial structures of the gauge theory. In particular the mod- 
uli are additive with respect to adding extra matter and increasing 
the number of colors. 
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1 Introduction 

Very recently there appeared new powerful methods of extracting effective 
superpotentials for a wide class oi M = 1 gauge theories. The proposal by 
Dijkgraaf and Vafa [1112], building upon earlier string theoretical construc- 
tions O H], links these superpotentials with quantities in random matrix 
models. The proposal have since been proven E] ■ The proposal has been 
extended to theories with fundamental matter ^ and there has been signifi- 
cant work on studying some features of the link with Seiberg-Witten curves 

lainiiiniiniiiaiiaiiii- 
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These curves made their appearance in the ground-breaking work of 
Seiberg and Witten [1^1 US] on the study of = 2 supersymmetric gauge 
theories. It turned out that one can describe the low energy dynamics of the 
gauge theory in terms of geometrical properties of the Seiberg- Witten curves. 
Subsequently it was reahzed jTHl US] that one could also study, within the 
same framework, deformations to A/" = 1 theories by adding a tree level su- 
perpotential. In this context one is led to the points (submanifolds) in moduli 
space where monopoles become massless and condense. At these points the 
Seiberg- Witten curve factorizes - it has only two single zeroes (branch points) 
in the case of complete factorization, where all monopoles condense. Once 
the form of the Seiberg- Witten curve at the factorization point is known one 
can calculate effective potentials for the deformed theory [ini UHl H] • 

The possibility of integrating in the glueball field S was realized in jH] 
and used to check the matrix model result following from the Dijkgraaf-Vafa 
proposal for deformed pure J\f = 2 theories (without fundamental matter). 
There the explicit factorization of curves without matter of ^H] was known. 
Factorization properties were also used in the context of SO{Nc) theories [20] 
and multi-trace operators fIT\ to obtain the appropriate effective superpo- 
tentials Weff{S). 

However once one adds fundamental matter to the theory, the factor- 
ization problem becomes exceedingly complicated (even in the case of U{2) 
with 1 flavour) and no explicit solutions are known. The aim of this paper is 
to use the Dijkgraaf-Vafa proposal linking superpotentials to matrix models 
and derive, from the random matrix model solution, explicit factorization of 
Seiberg- Witten curves of U{Nc) theory with Nf < Nc fundamental matter 
fields. 

If it were not for the Dijkgraaf-Vafa correspondence, we would not expect 
that any analytical solutions to this problem could exist. The fact that they 
can be obtained in this way shows that the random matrix model with matter 
is intimately linked to fine details of the geometry of the appropriate Seiberg- 
Witten curves. 

In addition we find that the matrix theory variables capture a surprising 
robust structure of the factorized Seiberg- Witten curve. The 'pure' Af = 2 
solutions and the new contributions of each flavor appear additively. All 
nonlinearity is concentrated in a single relation involving the scale of the 
gauge theory A and the matrix theory variables. Moreover we find that 
the 'integrating-out' equations of pure M = 2 theories reappear, in terms of 
matrix variables, in the context of A/" = 2 theories with fundamental flavors. 
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This arises naturally in the matrix model but is rather unexpected from a 
gauge theory perspective. 

In section 7 we collect the final results of the paper. 

The plan of this paper is as follows. In section 2 we describe the field 
theoretical ingredients in more detail. In section 3 we introduce the setup for 
the Dijkgraaf-Vafa proposal with matter fields. We then go to use the orthog- 
onal polynomial method to rederive directly the solution of the factorization 
for pure M = 2 theory, which is an ingredient of the expressions for theories 
with matter. In sections 5 and 6 we derive the linearity in couplings of the 
superpotential and the expressions for factorization. Section 7 contains the 
main result of the paper. In section 8 we give some specific examples and 
close the paper with a discussion. Three appendices contain some technical 
details. 



2 Field theory considerations 

It has been known for a long time that M = 2 supersymmetric U (N^) gauge 
theories with Nf < Nc massive fiavors has a Coulomb branch that is not lifted 
by quantum corrections This quantum moduli space is Ai'c-dimensional, 
parametrized for example by Uk = (||tr<I>'^^ with k < Nc- At each point of 
the moduli space, the low energy theory is described by an A/" = 2 effective 
abelian [7(1)^'= gauge theory. All the relevant quantum corrections in the 
IR can be recast in terms of the period matrix of a particular meromorphic 
one-form of the auxilliarly complex curve — the Seiberg-Witten curve ^3 
[TE\ [T71 ITS] , or more precisely a family of genus Nc — 1 hyperelliptic Riemann 
surfaces: 

Nf 

= Pr,Sx, u,Y - 4A2^-^/ + ^i)^ (1) 

i=l 

with 

PnA^^uu) = (det(x/-<l>)) = (2) 

0=0 

The coefficients Sa are polynomials of the m^'s parameterizing the Coulomb 
branch: 

a 

+ kSa-kUk = (3) 
fc=0 

So = 1 , uo = (4) 
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One can deform this M = 2 theory to a A/" = 1 gauge theory by adding a 
tree level superpotential: 

Wtree= E 9p ' ■ (5) 

p=l P 

The classical vacuum structure is given by all possible distributions of the 
Nc eigenvalues (f)k of $ amongst the critical points aj of the potential. 
This corresponds, at the classical level, to a breaking of the U{Nc) gauge 
symmetry to the gauge group ll^Ji '"'U{Ni) with J2i^=i^ = where the 
iVj's are the nonzero multiplicities of the eigenvalues. For the purposes of 
this paper we can safely set Qn^+i = as we will be considering the case with 
no breaking of gauge symmetry (complete factorization). 

Turning to the quantum picture, the presence of this superpotential will 
lift the quantum moduli space, characteristic of the M = 2 Coulomb phase, 
except for the codimension n submanifolds, where n mutually local mag- 
netic monopoles become massless. These are the M = 1 vacua solving the 
F-flatness and D-flatness conditions and are characterized by a monopole 
condensate of the massless monopoles. This is believed to produce, by the 
dual Meissner effect, the expected confinement of the electric M = 1 theory. 
Hence, the final quantum theory is described at low energies by a A/" = 1 
f/(l)^^-" gauge theory. These f/(l)'s can be thought of as the U{V) C U {Ni) 
of the classical theory. The A/" = 1 SU (Ni) part of the theory confines, has 
a mass gap and is characterized by a gaugino condensate. 

These Af = 1 vacua, being the codimensional n submanifolds of the M = 2 
Coulomb branch, where n mutually local monopoles become massless, are 
parameterized by the sets of moduli {u{^'^^'} where the Seiberg-Witten curve 
factorizes, i.e. the r.h.s. of ((T]) has n double roots^ and 2{Nc — n) single roots: 

i=l 

Moreover it is shown in that the reduced curve, 

= F2(N,-n) (7) 

captures the full quantum dynamics of the TV = 1 U{1)^''~"' low energy 
theory. 

^For simplicity we are assuming that higher order roots are not occuring. 
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The effective superpotential of the theory deformed by ©, where we set 
gNc+i = 0) is obtained by plugging the solutions parametrized by Nc—n 

parameters into the tree level superpotential: 

Nc 

'^eSi = Y.Sp<'-^ (8) 
p=l 

Then (jS} should be minimized with respect to the N(. — n parameters. 



Complete factorization 

In this paper we will be interested in the case where N^. — ! mutually local 
monopoles condense. This corresponds to a complete factorization of the 
Seiberg-Witten curve: the vacua form a 1 dimensional submanifold such 
that the curve has only 2 single roots and Nc — 1 double roots: 

PNXx,ut'-?-^^^'''~''^X{{x + m,) = ix-a)ix-b)Hl_,ix) (9) 

i=l 

The main goal of this paper is to find explicit expressions for the moduli 
ul"''^^' where complete factorization occurs. Let us briefly review the solution 
for the pure M = 2 Yang-Mills case found by Douglas and Shenker jTH] some 
time ago using Chebyshev polynomials. 

Their solution factorizing the curve Piq^{x^Uk)'^ — Ah?'^" is^ 




Note that the one dimensional submanifold where the U{Nc) curve com- 
pletely factorized, is parametrized by ui = tr$. These results can be easily 
restricted to the SU{Nc) case, by putting exphcitly Ui = 0. All parameters 
are then uniquely fixed in terms of the only scale of the theory A. 

Integrating in S 

The quantum JV = 1 effective potential generated by the tree level potential 
is obtained by minimizing (jSJ 

W,jf{A,u,) = j:9p<'H^,u,) (11) 

p=i 

^Adapted to the U{Nc) gauge group P]. 
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with respect to Ui. 

Here we follow an alternative route taken in jH], and integrate in S by 
performing a Legendre transformation with respect to log A^^'^~^-f . The su- 
perpotential is then given by: 

Weff{S,ui,Q,A) = S log A^^'^-^f + Weff{S,ui,Q) = 

= ^logA^^-^/ - 51ogfi2^-^/ + E^7p«r*-(fi,«iXl2) 

p=i 

Note that the only A-dependence is in the linear Slog A^^'^"^-^ term. Inte- 
grating out S forces Q = A and brings us back to (fTT|) . In order to get the 
effective potential iye//(A, S) one has to integrate out both Q and Ui. 

Our strategy for identifying the factorization parameters m^"^*' for the 
theory with matter is to rewrite the random matrix expression in the form 
given by the last two terms of (fT^ and to read off the appropriate gauge 
theoretic moduli m^"^*' . The first term in (jl2|) could also be absorbed into the 
matrix model expressions if appropriate rescalings of the integration measure 
of the matrix model was made. We will not do this here. 

For later reference let us quote the equations of motion for the pure M = 2 
gauge theory: 



These follow directly from ()12p and (fTUI) . 

3 The Dijkgraaf-Vafa proposal with funda- 
mental matter 

According to the Dijkgraaf-Vafa prescription the the perturbative part of the 
superpotential can be expressed as [1, 

WefAS) = N f^^=f^ + (15) 
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where the are defined through the matrix integral [7| 

e-|^-^x=2(5)-f^x=i(^^)+- = y /}$/)g./}Q.e-f (16) 

where we included the coupling of fundamental matter to the adjoint field in 
accordance with M = 1 supersymmetry. 

The matter fields in the matrix model ()16p appear only quadratically and 
hence may be integrated out giving 

1 det (m + $)-' e-f = Z ■ (det (m + $)"') (17) 

where Z is the partition function of the matrix model without matter. For 
the complex matrix model it is well know that Z will contribute only to T^=2 
and not to ^x=i- Using large N factorization we also have 

(det(m + ^Y^) = ^e-'°s'i^*('"+*)) = (e^^^ iog{m+#)^ ^ ^-(tr Iog(n^+^.)) ^^g^ 

We see that the matter determinant appears at subleading order in w.r.t 
the tree level potential. This has the important consequence that the saddle 
point solution (eigenvalue density p(A)) and hence J-'x=2 will not be influenced 
by the presence of matter. The x = 1 contribution is then given by 

Nf 

^x=i = E / ^V(A) log(m, + A) (19) 
1=1 

In the following section we will evaluate the J-'x=2 piece of the partition 
function Z appearing in ()17|) using the method of orthogonal polynomials. 
As a byproduct we will rederive the factorization expression for the pure 
gauge theory case. Then in section 5 we will start investigating the matter 
contribution (jl9|) . 

4 Orthogonal polynomials and J-x=2 

In this section we will use the method of the orthogonal polynomials to study 
thoroughly the one cut solution. A very nice introduction to this powerful 
method can be found in [221 ■ We will show that all results obtained from 
the field theory analysis appear naturally in this setting from random matrix 
computations (compare 0). 
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The orthogonal polynomials associated to the matrix model with poten- 
tial ^tr V^($) satisfy the recursion relation 

sPnis) = P„+i(s) + T„P„(S) + RnPnMs) (20) 

In the large N limit the recursion coefficients T„, Rn can be taken to be 
continous functions of u = S ■ {n/N) ^ S ■ x. In addition, the equations of 
motion of the matrix model with general potential Wtreei^) become purely 
algebraic [211221: 

^Lv'(, + ^ + T(S.)] = S. (21) 

<fa 1 ,„ ; ^ \ = (22) 

2m z \ z J 

When X = 1 we will denote R{S) by R and T{S) by T. These two variables 
are related to the endpoints a and b of the support of the matrix eigenvalue 
distribution through 

^ a + b , . 

T = ^ (23) 

(a -bf 

R = (24) 
For a general matrix potential V($) = Y^Qp^^^, one obtains easily: 

u ^ Sx = Y.9p E l[2q){ ) RiSxmSxy-"^ (25) 

V - j:9,E^(l](^-!]Risxmsxr-''^-' = o, (26) 



p>l q=0 



P X'^Q J \(1 J 



defining u and v for later convenience^. In terms of these polynomials, the 
matrix model partition function takes a very elegant form: 

Z = Ar!/i^e'^'/o''^-(i-)i-«(^-). (27) 

where ho is the integral 

ho = r dse-^^^'^ (28) 



^ Apart from the present section u and v are always taken at a: = 1. 
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From equation (j27j) . one can easily extract the x = 2 contribution (we will 
comment on the piece below) 

jr^^2 = -S^ C dx{l - x) In R{Sx) (29) 





/o 

Performing an integration by parts leads to: 

-S du In R{u) + duu In R{u) (30) 

JO JO 

Following the Dijkgraaf-Vafa prescription, the contribution to the gauge the- 
ory effective potential is proportional to: 

!^<4iiidfi (31) 



dS Jr{q)=o R 

The remaining integral is rather tricky. The integrand is a function of two 
variables R and T, tied together by the highly nonlinear constraint v = 0. 
This complicates the integral, and moreover spoils the manifest linearity in 
the couplings Qp, appearing in the matrix potential. It is convenient, however, 
to treat the one-dimensional integral as an integral of a 1-form over the path 
{v = 0} in the two dimensional R-T plane. Then one can rewrite the integral 
as an integral over a closed 1-form u, such that the original integral remains 
unchanged. One can then deform the contour v = into a piece with R = 0, 
while integrating T from T(0) to T{S) and another piece, keeping T fixed at 
T{S), and integrating R from to R{S) (see fig. 1): 

^ ' dR= ^ ' ' dR = uj= UJ+ uj 



'r{0)=0 R Jv=0 R Jv=0 JT=const J R=const 

(32) 

An extension of the integrand of (j^T|) to a closed one-form can be found to 
be 

u 

uj = -dR + vdT. (33) 
R 

Note that the extra piece doesn't give a contribution to the original integral, 
which is taken over v = 0. 

Performing the integral is now straightforward. The first piece, integrat- 
ing over T and constraining i? to gives: 



j^^^^ vdT = J2^{TP{S)-TP{0)). (34) 
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TA 



V v=0 



R 

Figure 1: The dotted line represents the original contour of integration, while 
the solid lines show the deformed contour in the T-R plane. 



The evaluation of the integral at T{0), can be cancelled (at the saddle point) 
against the Hq , appearing in the partition fuction (|77j) . For the second piece, 
one fixes T at T{S), and integrates over R: 



R{S) u 

-dR 
R 




j^qrpp-2q 



(35) 



From now on until the end of the paper we will denote by R and T the 
recursion coefficients defined at S. Bringing the two contributions together, 
leads to the following result: 



Weff = -N,S\nR 



p>l P q=0 




j^qrpp-2q 



(36) 



This is exactly, term by term, the field theory result for pure M = 2 U (iVc) 
gauge theory ((O and (fTUjl ). provided we identify field theory variables and 
matrix theory variables in the following way: 



R = n^ 



T 



V2_ 



(37) 



The first identification is extracted from the linear term in S, the second part 
follows from the term linear in gi. As it should, the equations of motion for 
R and T (j25|) - (j26|) . are mapped, under this identification, to the equations 
obtained from integrating out Q and Ui ()13j) - (fT^ . 
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From (j36|) we may read off the coefficients of Qp. 

U^^R, T) = i '2 ( ) ( ) (38) 

Under the identification ()37|) and setting Q = A these coincide exactly with 
the factorization solution (fTUj) of Douglas and Shenker for pure Af = 2 gauge 
theory. For theories with matter, (j38p will be just a part of the final result, 
and an identification between R, T and gauge theoretic quantities will have 
to be made only after the J-'x=i contribution is evaluated in the following 
sections. 



5 The J-x=i matter contribution 

In the previous section we have rederived the result that the first piece in 
()15|) can be recast in the form 



" dS 



-S\ogR + Y.9pU;'''\RM (39) 



and the equations of motion for R and T derived from ()39|) are exactly the 
random matrix saddle point equations 

The first of these can be interpreted as the formula for integrating in S. The 
linearity of (jH^ll in the random matrix couplings is essential for identifying the 
Ijpure^j^^ rp^ with the point in AT = 2 moduli space where the Seiberg-Witten 
curve y"^ = P^X^i'^k)'^ ~ 4A^^= factorizes. 

In the following section we will recast the random matrix expression ()19|] 
for J^x=i in the same way: 

Nf 



^ r n 

•^x=i = E [^logL(/?,T,m,) +5:<7,W7«^^^(/2,T,m,)] (42) 

i=l 
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Let us comment on an ambiguity, leading to a unique determination of a 
crucial term to obtain the correct final result. Once we start interpret this ex- 
pression for T^=\, as a part of the effective gauge theory potential, we should 
be able to integrate out R and T. The equations of motion thus obtained, 
should be consistent with the random matrix saddle point equations. But, 
as explained in section 3, adding fundamental matter in the gauge theory 
does not modify the saddle point equations of the matrix model. So it is 
necessary to impose that the equations of motion for R and T, derived from 
the effective potential, with the T^=\ contribution are consistent with ()40p 
and ()41|) . This fact is highly non-trivial and reveals an unexpected relation 
between the factorization of Seiberg-Witten curves with and without matter. 

We will show that the addition of T-^=\ to the superpotential does not 
change the equations of motion for R and T provided we add an extra term 
proportional to v. From the matrix model perspective nothing changes, while 
the extra term turns out to be crucial to obtain the correct gauge theoretical 
interpretation of our results. We stress that this does not alter the Dijkgraaf- 
Vafa prescription. It merely solves an ambiguity that arises when interpreting 
the matrix model variables directly in terms of gauge theoretical quantities. 

Let us now focus on a single summand of ()19p. In this case it seems to 
be difficult to use the orthogonal polynomial method, which was the most 
straightforward way of deriving the T^=2 contribution. Here it is more con- 
venient to use the saddle point expression for the eigenvalue density of the 
single-cut solution: 



where M{x) is a polynomial which is expressible in terms of the random 
matrix potential through 



The features of interest of the above expression are (i) it is linear in the 
explicit dependence on the couplings (ii) the coefficients of are universal 
functions of the endpoints a, h and hence of the variables R and T, (iii) for 
given (yfp, its coefficient is a polynomial in x of order p — 2. 

In order to obtain an explicit dependence on R and T let us perform the 




(43) 




(44) 
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change of variables 

x = ^{a + b) + ^^^ij = T + 2\^iP (45) 
Then the contribution to J-'x=i of a single flavour is given by 

R-^J\ dip^Jl - ■ M(V') ■ log (m + T + 2Vr^p) + vf{R, T) (46) 

As noted before, one should take into account a possible addition of v mul- 
tiplied by any function /(-R, T) (since v = is an equation of motion of the 
random matrix model). 



6 Factorization formulas 

We will first derive the formulas involving gi and in particular this will 
allow us to fix uniquely the function f{R, T). Also this will make the general 
structure more transparent. Then we will derive the results for arbitrary gp. 



Contribution of gi, g2 

To this order M{iIj) = g2 and the formula 



gives 



g2R 



log(m + T) + ^ j^^ d'lp^l- ^/'2iog M + 2 



R 



m + T 



(47) 



This can be evaluated to give 
92R log 



m + T + J{m + Ty -AR 



+ (m + T)- 



m + T - J{m + T)^ -AR 1 



AR 



(48) 

At this stage we should identify the coefficient of the logarithm with 5* (here 
this is trivial since to this order the equation of motion for R is just S = g2R, 
but later we will see that this property will hold in general). Thus we are 
left with 



5 log 



m + T + J{m + Ty - AR 



+92 



m + T 



R 



^ ,m + T - ^{m + ry -AR) ^ 

(49) 
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which indeed has the form of (j42j) . Now we require that the saddle point 
equations remain consistent with the T^=\ equations of motion. This detem- 
ines uniquely the term 

v{g^, g2, R, T) ■ f{R, T) = {g^ + g^T) ■ f{R, T) (50) 

Indeed the requirement that integrating out T from the sum of and (jMIjl 
gives V = gi + g2T = fixes f{R, T) uniquely to be 

/(i?, T) = -]^[m + T- ^{m + TY-AR) (51) 

This extra term does not change the equation for R, which is consistent with 
the saddle point equations. This function will stay unchanged in the general 
case. We may now read off the final expressions for the mass dependent 
contributions to Mi and U2: 

Ur''"{R, T,m) = -i (^m + T - ^(m + T)^ - AR^ (52) 
Ur'''''{R,T,m) = !^(m + T-^(m + T)2-4i?)-| + 

~(^m + T- ^{m + T)2 - 4i?j (53) 



Arbitrary gp 

We will now extend the previous considerations to the calculation of arbitrary 
■jjmatter ^ The general structure will remain unchanged. The function f{R,T) 
multiplying v will remain unmodified (as it should). Also the coefficient of 
the logarithm will turn out to be exactly S. 

In appendix A we derive the following expression for the polynomial 
M(V'): 

M{^P) = Y^gpJ2cp^r.r (54) 

p>2 n=0 

where 

r p-n-2 ] 

l 2 J /o;. \ / „ 1 \ 

j^kr^p-n-2-2k ^55^ 
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So the integral (j46j) can be rewritten as 

p-2 , 2 



p>2 n=0 

2 '•I 



+ - J ^ diP^Jl - ^2 . . log M + 2 



m + T 



(56) 



We now have to distinguish two cases. 

n odd: Then the first integral vanishes and we will denote the second 
integral by fn{z). As discussed in appendix B, fn{z) is essentially a poly- 
nomial in z of order (n + l)/2 (divided by {z — l)"/2+i^ when expressed in 
terms of the variable 

m + T 



2R 



m + T + J{m + Ty - AR 



(57) 



Appendix B contains a general formula for fn{z)- Explicit expressions for 
some specific cases are shown in table 1 in section 7. 

n even: In this case the first integral is nonvanishing. Moreover the 
second integral involves a logarithm with the same coefficient as log(m + T) 
in the first integral. Together they combine to give 

. - ^ 'm + T + J(m + T)2 -4i?' 



'21 



p>2 



1=0 



l+l 



21 



} ■ lo£ 



(58) 



It is shown in appendix C that the coefficient in curly braces is exactly equal 
to S. The second integral with the logarithmic part subtracted out has again 
a simple polynomial structure (see appendix B for details). 
At this stage we arrive to the analogue of (jl^ : 



5* log 



m + T + J{m + Ty -AR 





V2 




^ ^ Cp^nRfni^^) 


/ P>2 


n=0 



Again we have to add to this the correction term 



v.f{R,T)^-Y^g,Vp-(m + T-^im + Ty 
p>i 



AR 



(59) 



(60) 



We checked explicitly for some cases that this term together with ()59|) gives 
equations of motion consistent with the random matrix constraints. 

The sum of (l59|l and (j60|l is now of the expected form (jl2|l . thus defining 
(i?, T, mi). 



^matter 
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7 Final results 

Putting all results together and putting in the term S" log A^^^"^^ 

as required from the Dijkgraaf-Vafa prescription gives a prediction from the 
matrix model side for the quantum effective gauge potential: 



Weff{S,T,R,A) = S\og ^ 



p>i 



i=l 



■ (61) 



This expression should be compared with the potential Weff{S, ui, fl, A) (see 
()12|l ). obtained from the field theory analysis. The relation between the 
parameters of the matrix model and the field theory is highly non-linear: 



u, = N,T--J2{mi + T- v'K + Ty - AR) (62) 

i=l 

^^^^-"^^ = vvTTT „ . (63) 



n:L'i l[mi + T + y'(m, + T)2 - 4R 

In order to obtain the final factorization formulae we integrate out S from 
We//(5,T,i?,A) : 

^2N^-N, ^ ^2N.~N, ^ (64) 

Combined with ()63|) this gives an expression for A in terms of R and T. The 
remaining part of (|6ip should then be compared with field theory result: 

Weff{A,m) = ^ <7,«f^*-(A,«i), (65) 
p>i 

where the uf/"'-, is the one parameter solution, factorizing completely the 
Seiberg-Witten curve for f/(iVc) SYM with Nf flavours {Nf < N^): 

= PnSx, Uk)' - 4A2^-^/ n(^ + (66) 

1=1 



Comparing with the results from matrix models gives an expression for the 

Nf 

^f^act. ^ N^lfPure^j^^ T)+J2 ^^""''(i?, T, 771,) (67) 

i=l 



Ufact.,^ 
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2(2-1) ^IV^; — 6(2-1)3/2 

f — 1 f (■y\ — 302^-652+32 

-3z^+9z-5 f I ^\ — 5252^-161022+15822-512 



16(2-1)2 J'iy^) 120(2-1)5/2 

J- / N _ -322+92-5 f _ 5252^-16102^ 

J4l,^; — 96(2-1)3 J5l^; — 3360(2-1)7/2 

f / _ -4823+16822-1762+59 f ^ _ 44102^^-1764023+2595622-168572+4096 
J6l^j — 1536(2-1)4 i7l^J — 40320(2-1)9/2 

Table 1: Examples of the functions fn{z) for small n. 



in terms of two parameters R and T, tied together with the constraint: 

DiVc 

nS I (m, + T + + T? - AR) 

For completeness, we recall the formulas 

ypure T) = ^ ^ ) ( ) (69) 



Ur''"{R, T,m) = --(m + T - ^(m + T)2 - AR ) (70) 



U^^t'iR, T,m) = J2 Cp,nRfn{z) -Vp-[m + T- ^{m + Tf - ARpl) 



n=0 



In the above formula the coefficients Cp^„ are defined in (|55|). while Vp is just 
the coefficient of Qp in the constraint v = (see eq. ()41|) ). Finally the func- 
tions fn{z) are computed in appendix B and depend on the variable z{R, T), 
given by (j37|) . In table 1 we present the explicit forms of the functions fn{z) 
for n < 7. 

In the final result (jU7|) we see that both the and Nf dependence is 
very simple. Moreover the contribution of each extra fiavor enters additively 
the expression for the moduli. Another curious feature is the appearance of 
the original factorization solutions for the pure M = 2 U (Nc) gauge theory. 
Increasing the number of colors does not change the expressions for ]J^°-*^<^'^ 
in terms of R and T. However the only nontrivial change is encoded in the 
expression for A in terms of R and T. 

If we wanted instead to obtain the effective potential W{S,K), we would 
have to integrate out R and T from ()61|) . On the field theory side it is 
very cumbersome how the structure of the Seiberg-Witten curve appears in 
the equations of motion for Ui and VL. On the matrix model, on the other 
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hand, the equations of motion for R and T appear naturally to be the same 
as ones obtained in the case without flavours. It seems that the particular 
combinations of Ui and Q, embodied in R and T, captures some nontrivial 
structure of the Seiberg-Witten curves. 



8 Some examples 

In this section we will study some examples to verify that the i*^"'^*" we 
obtained from random matrix models, do factorize the appropriate Seiberg- 
Witten curves with fundamental matter. 



U{2) with 1 flavour 

In this case we have 



/act. ^2T-^(m + T ~ ^{m + Tf - AR^ (72) 



"2 



9 R2 

^3 ^ (74) 

m + T + J{m + T)2 - 4i? 



We have verified that with the above choices, the discriminant of the Seiberg- 
Witten curve 

= (^"^ - uix -U2 + ^uj^ - 4A^(x + m) (75) 

vanishes identically, which in the special case of 2 colours proves factorization. 
Note that for general iVc, complete factorization is a much stronger condition 
than the vanishing of the discriminant. 



SU{Nc) with 1 flavour 

In order to consider SU{Nc) theory we have to impose the constraint that 
Ui — 0. Then the parameter T can be expressed in terms of R and the mass 
m of the additional flavour. Namely we have 

ui = NJ' -]^[m + T- yj{m + Tf - 4R^ = (76) 
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which gives 



T = ^ (77) 

2(iV,-l) 

Now R is hnked directly to the scale of the gauge theory through 

A^^-^ = (78) 

{m + T + ^{m + T)2 - ARj /2 

where the expression (|77p should be used. The remaining formulas remain 
however quite complicated functions of R, m and N^. This is in marked 
contrast to the case of pure gauge theory without fundamental matter where 
the passage from U{Nc) to SU{Nc) is very simple, and enters linearly. 

It is interesting to look at the m —>■ oo limit. Then T ^ as expected for 
a pure A/" = 2 SU{Nc) theory, while ^ becomes A^^-^-i = R^-/m. Recall 
that in pure SU{Nc) theory R had the interpretation of A^^^^. Hence we 
obtained the correct field theoretic matching of scales 

A^^^-^m = AJ^^, (79) 

Moreover it is easy to check that then the W^'"**'^''(i?, T, m) give a vanishing 
contribution as the functions fn{z) when z ^ oo (see appendix B). The 
above behaviour is quite clear from the Seiberg-Witten curve perspective. It 
is reassuring that it could be also obtained in a simple way from the random 
matrix formulas. 



U{7) with 3 flavours 

As a final check of the formulas we considered U{7) theory with 3 flavours 
with masses mi = 1, m2 = 2, = 3. For the random choice of parameters 
T = and R = 0.2 we flnd A = 0.31643, while the polynomial Pjix) = 
x^+0.45018x''-1.3447a;^-0.5382a;^+0.5048x^+0.16048a;2-0.045a;-0.00704. 
For the curve 

y2 = p.^(a;)2 _ 4^11(3; + + 2)(x + 3) (80) 

we flnd single zeroes at x = ±0.8944 and a series of 6 double zeroes in 
between. 
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9 Discussion 



In this paper we used the Dijkgraaf-Vafa proposal hnking random matrix 
models with fundamental matter and superpotentials for obtaining explicit 
formulas for the complete factorization of Seiberg-Witten curves for U{Nc) 
theories with Nf < Nc flavours. These points in the moduli space, forming 
effectively a 1-parameter manifold, correspond to condensation of all species 
of monopoles. As a byproduct we obtained formulas for the solution of the 
random matrix model with matter with an arbitrary polynomial superpoten- 
tial. 

In order to identify the points in moduli space where the Seiberg-Witten 
curve factorizes we recast the random matrix solution in a way that exhibits 

(i) linearity in the couplings gp of the deforming tree level superpotential, 

(ii) the whole dependence on the glueball superfield S could be written as a 
linear coupling of S* to a logarithmic expression. The first property allowed 
us to identify the moduli space parameters of the factorized curve Up as 
the coefficients of Qp, while the second property is exactly the one found in 
'integrating-in' 5* and thus gave the expression for the gauge theoretic scale 
A in terms of random matrix model quantities. 

The fact that the above procedure works is yet another argument for the 
Intriligator-Leigh-Seiberg linearity principle [21] and the validity of 'integrating- 
in'. In addition it shows that the matrix model of the Dijkgraaf-Vafa pro- 
posal captures quite detailed properties of the field theoretical Seiberg-Witten 
curve. In fact we found it surprising that any analytical description of the 
very nonlinear complete factorization property could be found for the case 
with matter. 

A curious feature of the random matrix formulas is that the solution for 
the Up for the pure M = 2 theory appears linearly in the complete expression 
for the theory with fundamental matter fields. The full 'nonlinearity' is 
encoded in the formula for A in terms of random matrix parameters. It 
would be interesting to understand this structure from the field-theoretical 
point of view. In addition the random matrix constraints expressed in terms 
of R and T don't change when adding fundamental matter. They have 
precisely the form of equations of motion for the pure Af = 2 theory. But 
now the mapping between R and T and field theoretical A and ui becomes 
complicated. Nevertheless, the question why the pure Af = 2 equations 
still arise in a disguised form for theories with fundamental matter poses an 
interesting question from the gauge theory perspective. 
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There are numerous issues that one could investigate further. The fac- 
torization properties of the pure Af = 2 curve are hnked with the concept 
of master field. This has been investigated in the context of associated ran- 
dom matrix theory in |25j . It would be interesting to study the factorization 
formulas obtained in this paper from a similar point of view, albeit it will 
surely be much more involved. 

Another interesting question would be to explore the mathematical struc- 
ture linking the random matrix model with matter with factorization proper- 
ties of the associated curves. In this paper we relied heavily on recasting the 
random matrix expressions guided by field theoretic ingredients such as the 
ILS principle and integrating-in, for which there is no real direct proof. It 
would be very interesting to uncover the mathematical interrelation between 
such seemingly unconnected topics as the factorization of SW curves and 
random matrix models. 

Finally we hope that the above results could be used for a more detailed 
investigation of the physics of U{N^ theories with Nf < flavours along 
the lines of [13120]. 

Acknowledgments RJ was supported by the EU network on "Discrete Ran- 
dom Geometry" and KBN grant 2P03B09622. YD was supported by an 
EC Marie Curie Training site Fellowship at Nordita, under contract number 
HPMT-CT-2000-00010. 



Since we use explicitly the form of the eigenvalue density in the variable ip, 



A Formula for M^ijj) 



let us perform the changes of variables x 
the definition (jSl) of M(x): 



T + 2^/R^|J, w = T + 2VR(f) in 



M{ip) 



1 



Res0: 



1 V'{T + 2^/R(f)) 



(81) 



2VR 



i=oo 



Using the power series expansion of the square root 




1 



1 



,2k 



(82) 
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it is straightforward to obtain the Laurent expansion of the function in (|8ip: 



OO OO 1 p— 1 



E^.EEr-5^a.E r / M (2v^)"t--V (83) 



/,2fc+n+2 « ^ \ J 
p>2 n=Ok=0 "r 1=0 \ / 

From this expression we may isolate the coefficient of 1/0 giving the result 
quoted in the text: 

= T.9pT.^P,nr (84) 

p>2 n=0 

\ p — n — 2 ] 



fc=0 



B Logarithmic integrals fn{z 

Here we will derive the explicit form of the functions fn{z) related to the 
logarithmic integrals 

4 = - d^yiT^.^«.log(i + 2a;^) (86) 

where x = \/R/{m + T). In fact the results simplify significantly if one 
reexpresses everything in terms of the variable z 

z = ^^^^ {m + T+ ^(m + T)2 - 4i?) (87) 

{x is expressed in terms of z as \/z — Ij z). We have to distinguish two cases: 

n odd 

The integral (jHUj) can be performed using a series expansion of the logarithm, 
integrating it term by term using 

and resumming. The result is 

A{z-l) 



5 + n 



;z2 



(89) 
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For odd n, fn{z) is expressed through elementary functions (see examples in 
section 7). It has the form of a polynomial in z of order (n + l)/2 divided 
by {z-l)^+\ 



n even 



The integral ()86|) can be again obtained using a resummation procedure. The 
result is 



-- ^^T{z + f I 1.1,— 2.3 + 



4(2 



^2 



(90) 



In this case the integral ()86|) involves a logarithm, which together with 
the log(m + T) forms the logarithmic function multiplying S in ()59|) . Thus 
to define the function /„(z) for even n we have to subtract from J^^^" this 
logarithm: 

/r"(^) = /r" - ^ f n V log ("-^\ (91) 



Its general form turns out to be a polynomial of order n/2 divided by 

{z — 1) 2'+^. In table 1 in section 7, for completeness we present the explicit 

forms of the functions fn{z) for n < 7. 



C Coefficient of the logarithmic term in ( 1581) 

In this appendix, we identify the coefficient of the logarithmic piece of the 
X = 1 contribution with S, as expected from field theory. From (j46p one can 
easily read off the coefficient of log(m + T): 

4i? £ ^yr^M(^) (92) 

Using the elementary integral (j88j) one obtains the expression: 

/9/\ 2-2'-2 
^REdpEi ]^P,2ij-Y (93) 

p>2 21=0 \ W ' + ^ 

Inserting the explicit expression for the coefficient Cp^2i, leads to: 

E * E ± f 2, : 1 ^fi«-T— (94) 



p>2 1=0 k=0 
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Changing to a new summation variable m = k + I + 1 gives the expression: 

E*lE^f,!iU?V^::"-,1^"r-" (05) 



This is exactly the expression for S, provided that: 

1 (2l\ (2m-2l-2\ (2m\ 

which can be verified. 
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